Single-cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10× Genomics Chromium single-cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human monocyte-derived dendritic cells infected with ZIKV at the single-cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN-dependent and -independent genes (the antiviral module). We modeled the ZIKV-specific antiviral state at the protein level, leveraging experimentally derived protein interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per-cell basis with experimental protein interaction data. IMPORTANCE: Zika virus (ZIKV) remains a public health threat given its potential for re-emergence and the detrimental fetal outcomes associated with infection during pregnancy. Understanding the dynamics between ZIKV and its host is critical to understanding ZIKV pathogenesis. Through ZIKV-inclusive single-cell RNA sequencing (scRNA-seq), we demonstrate on the single-cell level the dynamic interplay between ZIKV and the host: the transcriptional program that restricts viral infection and ZIKV-mediated inhibition of that response. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool for gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Frequency-potency analysis of IgG+ memory B cells delineates neutralizing antibody responses at single-cell resolution.

Identifying individual functional B cell receptors (BCRs) is common, but two-dimensional analysis of B cell frequency versus BCR potency would delineate both quantity and quality of antigen-specific memory B cells. We efficiently determine quantitative BCR neutralizing activities using a single-cell-derived antibody supernatant analysis (SCAN) workflow and develop a frequency-potency algorithm to estimate B cell frequencies at various neutralizing activity or binding affinity cutoffs. In an HIV-1 fusion peptide (FP) immunization study, frequency-potency curves elucidate the quantity and quality of FP-specific immunoglobulin G (IgG)+ memory B cells for different animals, time points, and antibody lineages at single-cell resolution. The BCR neutralizing activities are mainly determined by their affinities to soluble envelope trimer. Frequency analysis definitively demonstrates dominant neutralizing antibody lineages. These findings establish SCAN and frequency-potency analyses as promising approaches for general B cell analysis and monoclonal antibody (mAb) discovery. They also provide specific rationales for HIV-1 FP-directed vaccine optimization.

Single cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells.

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10x Genomics Chromium single cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human moDCs infected with ZIKV at the single cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN dependent and independent genes (antiviral module). We modeled the ZIKV specific antiviral state at the protein level leveraging experimentally derived protein-interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per cell basis with experimental protein interaction data. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool to gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Investigation of COVID-19 Outbreak among Wildland Firefighters during Wildfire Response, Colorado, USA, 2020.

A COVID-19 outbreak occurred among Cameron Peak Fire responders in Colorado, USA, during August 2020-January 2021. The Cameron Peak Fire was the largest recorded wildfire in Colorado history, lasting August-December 2020. At least 6,123 responders were involved, including 1,260 firefighters in 63 crews who mobilized to the fire camps. A total of 79 COVID-19 cases were identified among responders, and 273 close contacts were quarantined. State and local public health investigated the outbreak and coordinated with wildfire management teams to prevent disease spread. We performed whole-genome sequencing and applied social network analysis to visualize clusters and transmission dynamics. Phylogenetic analysis identified 8 lineages among sequenced specimens, implying multiple introductions. Social network analysis identified spread between and within crews. Strategies such as implementing symptom screening and testing of arriving responders, educating responders about overlapping symptoms of smoke inhalation and COVID-19, improving physical distancing of crews, and encouraging vaccinations are recommended.

Mumps Outbreak in a Marshallese Community - Denver Metropolitan Area, Colorado, 2016-2017.

In January 2017, the Colorado Department of Public Health and Environment (CDPHE) identified four epidemiologically linked cases of mumps among persons from a Marshallese community who were members of the same church in the Denver metropolitan area. During 2016-2017, sizable outbreaks of mumps reported in Arkansas, Hawaii, and Washington also affected the Marshallese population (1). CDPHE, the Tri-County Health Department (TCHD), and Denver Public Health collaborated to conduct an outbreak investigation during January-March 2017 using active and passive surveillance that identified 17 confirmed and 30 probable cases. Public health actions included conducting measles-mumps-rubella (MMR) vaccination clinics at local Marshallese churches; these resulted in the vaccination of 126 persons with ≥1 doses of MMR vaccine. Implementation of active surveillance and support from local Marshallese church leaders in promoting vaccination programs likely contributed to interruption of the outbreak.

Loop Sequence Context Influences the Formation and Stability of the i-Motif for DNA Oligomers of Sequence (CCCXXX)4, where X = A and/or T, under Slightly Acidic Conditions.

The structure and stability of DNA is highly dependent upon the sequence context of the bases (A, G, C, and T) and the environment under which the DNA is prepared (e.g., buffer, temperature, pH, ionic strength). Understanding the factors that influence structure and stability of the i-motif conformation can lead to the design of DNA sequences with highly tunable properties. We have been investigating the influence of pH and temperature on the conformations and stabilities for all permutations of the DNA sequence (CCCXXX)4, where X = A and/or T, using spectroscopic approaches. All oligomers undergo transitions from single-stranded structures at pH 7.0 to i-motif conformations at pH 5.0 as evidenced by circular dichroism (CD) studies. These folded structures possess stacked C:CH(+) base pairs joined by loops of 5'-XXX-3'. Although the pH at the midpoint of the transition (pHmp) varies slightly with loop sequence, the linkage between pH and log K for the proton induced transition is highly loop sequence dependent. All oligomers also undergo the thermally induced i-motif to single-strand transition at pH 5.0 as the temperature is increased from 25 to 95 °C. The temperature at the midpoint of this transition (Tm) is also highly dependent on loop sequence context effects. For seven of eight possible permutations, the pH induced, and thermally induced transitions appear to be highly cooperative and two state. Analysis of the CD optical melting profiles via a van't Hoff approach reveals sequence-dependent thermodynamic parameters for the unfolding as well. Together, these data reveal that the i-motif conformation exhibits exquisite sensitivity to loop sequence context with respect to formation and stability.

Calorimetric and spectroscopic investigations of the binding of metallated porphyrins to G-quadruplex DNA.

BACKGROUND: The human telomere contains tandem repeat of (TTAGG) capable of forming a higher order DNA structure known as G-quadruplex. Porphyrin molecules such as TMPyP4 bind and stabilize G-quadruplex structure. METHODS: Isothermal titration calorimetry (ITC), circular dichroism (CD), and mass spectroscopy (ESI/MS), were used to investigate the interactions between TMPyP4 and the Co(III), Ni(II), Cu(II), and Zn(II) complexes of TMPyP4 (e.g. Co(III)-TMPyP4) and a model human telomere G-quadruplex (hTel22) at or near physiologic ionic strength ([Na(+)] or [K(+)]≈0.15M). RESULTS: The apo-TMPyP4, Ni(II)-TMPyP4, and Cu(II)-TMPyP4 all formed complexes having a saturation stoichiometry of 4:1, moles of ligand per mole of DNA. Binding of apo-TMPyP4, Ni(II)-TMPyP4, and Cu(II)-TMPyP4 is described by a "four-independent-sites model". The two highest-affinity sites exhibit a K in the range of 10(8) to 10(10)M(-1) with the two lower-affinity sites exhibiting a K in the range of 10(4) to 10(5)M(-1). Binding of Co(III)-TMPyP4, and Zn(II)-TMPyP4, is best described by a "two-independent-sites model" in which only the end-stacking binding mode is observed with a K in the range of 10(4) to 10(5)M(-1). CONCLUSIONS: In the case of apo-TMPyP4, Ni(II)-TMPyP4, and Cu(II)-TMPyP4, the thermodynamic signatures for the two binding modes are consistent with an "end stacking" mechanism for the higher affinity binding mode and an "intercalation" mechanism for the lower affinity binding mode. In the case of Co(III)-TMPyP4 and Zn(II)-TMPyP4, both the lower affinity for the "end-stacking" mode and the loss of the intercalative mode for forming the 2:1 complexes with hTel22 are attributed to the preferred metal coordination geometry and the presence of axial ligands. GENERAL SIGNIFICANCE: The preferred coordination geometry around the metal center strongly influences the energetics of the interactions between the metallated-TMPyP4 and the model human telomeric G-quadruplex.